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ABSTRACT 



Considering the absorption of terrestrial radiation by carbon dioxide 
and water vapor, coefficients of eddy conductivity are computed for heights 
in the range 4 to 320 meters from data collected during the Great Plains 
Field Program (August and September 1953.) Power-law profiles with height 

c, 

are fitted to the computed coefficients and values obtained from the fitted 
profiles correlated with a stability parameter. 
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Introduction , 



An important step forward in the theory of the turbulent transfer of 
properties in the lower layers of the atmosphere was taken vmen the con- 
cept of Au st.au schkoeffizienten or "exchange coefficients" was proposed 
by Wilhelm Schmidt in 1925. Since that time several attempts have been 
made by various workers to measure these coefficients through experi- 
mental observations in order to test the theory of turbulent exchange. 

In .1950 Cooley / 5_/ made determinations of eddy conductivity near the 
earth's surface from data gathered at Manor, Texas. About the same 
time Jehn and Gerhardt /9_/ made measurements of eddy heat flux up to 
about 120 meters at Manor and Gila Bend, Arizona. In addition to these 
were the observations of Cramer and Record /6 _/ near the earth's surface 
in 1353, and a small number of observations due to other workers in the 
period 1955-1957. These observations have been summarized by Priestly 
t and no attempt will be made here to review all of the results. 

\ 

The primary reason that the number of such observations has been 
small is that the accurate determination of micrometeorological parameters 
is exceedingly difficult and very expensive . For this reason the quality 
and volume of such data obtained in the past have been far below the 
level of desirability. 

The publication of the "Proceedings of the Great Plains Turbulence 
Field Program, 1 August to 8 September 1953, O'Neill, Nebraska" made 
new and better data available for the study of the turbulent transport of 
properties in the atmosphere's lowest layers. In order to obtain nev«” 
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I 





evidence on the nature of the turbulent heat and momentum transports 

and their relationships to one another, it was originally decided to 

initiate thesis projects at the United States Naval Postgraduate School. 

These projects to be undertaken were to determine the coefficients of 

(1) eddy conductivity, K„, (2) eddy viscosity, K.,, and (3) the ratio 

xl M 

of the tv/o coefficients, , from this new micrometeorological data. 

H M 

In particular this present study involves the determination of the coeffi- 
cient of eddy conductivity and the examination of its variation with 
height and stability from the values determined. 
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The Great Plains data. 



The data collection program for the Great Plains project was brohen 
down into seven time-observation periods with a multitude of parameters 
being obser/ed at nearly equally spaced observation times vrithin each 
period . 

« 

For this study the seventh observation period (7 and 8 September) 
vvis used and was considered to be the most suitable of the periods on 
the basis of a combination of two factors; (1) completeness of data uti- 
lized in the study and (2) a minimum of advection during the period. As 
an aid in the interpretation of the results the following quotation / 
is provided concerning the synoptic conditions in the area during the 
seventh period: 

On 7 September, the high-pressure area east of Nebraska moved 
eastward v/hile a trough of low pressure approached from the west. The 
winds were southerly and initially light, becoming 15 to 20 mph. The 
skies were generally clear over Nebraska with some scattered cirrus 
during the morning of the 8th. Cumulonimbus clouds and lightning were 
visible to the north of the test site from the afternoon of the 7th, until 
the early morning of the 8th. The seventh and last general observation 
period began at 0405 CST on 7 September and continued until 1305 GST 
on B September under relatively favorable conditions. 

Of the many tables of data included in /l^/, the following summary 
shows those which were used in this study: 

T able number ^ 

i.Ia 7 September, 1900, CST, West and 

South components of the geostrophic 
wind at indicated pressure levels. 
Values estimated from local synoptic 
charts , 

^Table where numbers of the form 1 . 1 etc . , refer to tables in /lO/ 
•'W'here times are omitted , data were available for each sounding 



A 

Data use d ^ 
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*5 ***> 

Table Number*’ Data Used 



1.3 


Station pressure and reicu e .... iuJ:y 


2 . la 


Soil temperature at O.S cr'. belovv tne 
surface of the ground (taken as being 
representative of the ground tempera- 
ture . ) 


4.2 *• 


Mean air temperature from a.C',,ir.i.ted 
thermocouples (Mil, 15-rr;Afiute raean- 
at2.0, 4.0, 8.0 and 1 € . 0 metej abo/e 
the ground , U CLA , average of Lix 
readings taken every 3 sfuiiutes at 
various levels up to 8.0' meter;;); 
shielded thermocouples (JK, average 
of 20 readings at each level during a 
5-minute interval at varioub lev’eis up 
to 6.4 meters) and ef-ieldecl therrnisters 
(UT , 10-second mear»r at \'irio'v: . levels 
up to .8 meters.) 


5.2a 


Mean cross product ot v’ercical eddy 
velocity component and air tempera- 
ture multiplied by density and specific 
heat of air for 3.0, 6.0 and 12 „0 meters 
above the ground, valuer nussir.g for 
1235 and 1435. 


6.4 


Radiosonde pressure, temperature and 
relative humidity at significant lev els . 


6 . 5 


Air temperature and mixing ratio ar alti- 
tudes from aerograph (L-2U) data, i 
at start , top level and end of lUglst ic 
indicated , 


7 9 3a 


Sensible heat transfer to the air iron, 
the earth-air interface a.ccord.ng to 
theoretical models of Suomi unu Lettau. 



During the observation period, data were recorded bi-hoiirly for e<.en 

i'ours of average local time (CST plus 35 minutes) . Sufficient time and 

viata v/isre available to make K„ computations for six soundings, 0t.3S, ObSS, 

1335, 1235, 1435 and 1635 on 7 September. Since no night-time aerograph 

observations vrere made, only daytime determinations were possible. 

4 



2,1. Instrumental errors . 



Contributors to the Great Plains project have given the following 
/alues for instrumental errors inherent in the data; 

I 

1 . Micrometeoro logical absolute temperature errors (Table 4 . 2) 

a) 0.2°C (UT) 

«b) 0 . 1°C (MIT and UCLA) 
c) 0.0 5°C (JH) 

2 

2. The L-20 aerograph configuration temperature and relative 
humidity errors . 

a) Temperature error: 

(1) 1°C absolute , for any level . 

(2) Relative error between levels, 0.1°C. 

b) Relative humidity error: 

(1) to 10% absolute, 

(2) Relative error between levels, 3%. 

In 29 cases a comparison of MIT temperature values at 16,0 meters 
to the L-20 data at 17 meters gave a mean difference, MIT minus L-20 of 
- 0,10 j, 0.11°C. 

For a complete picture of the instrumentation, analysis of instru- 
rRsntal errors and a description of the test-site the reader is referred to 

ZW- 



1 

•^""For an interesting comparison of observed micrometeoro logical tempera- 
tures see fig , 1 , 

L~20 configuration utilizes the KS-4 Kollsman aerograph, which is 
an automatic device for the automatic recording of temperature and rela- 
tive humidity , 




temperature observations for the 0635; sounding. 



o 




Background, methods and equations. 



The coefficient of eddy conductivity, K^j., is defined by; 



Kh - - 



Q 



pCp c © 



v/here Q is the turbulent heat transfer in the vertical. It may be observed 
from this definition that K„ is always a positive quantity. Wfien ^ is 
negative (unstable conditions) , Q will be positive and directed upward. 
Conversely for stable conditions >o) > Q will be negative and directed 
downward . 



In this study starting values of Kj^ at 4 and 8 meters are coropated 
by using eq. (1) with Q given by the average value in the layer 3-12 
meters. For values above the starting level the following method is utilized, 
The horizontal turbulent exchange of heat is usually very small 
compared to that in the vertical, and vertical advection is usually con- 
sidered to be negligible in comparison to horizontal advection. The fact 
that clouds did not occur in the layer under study rales out a latent heat 
contribution , and the small height -variation of density can be overlooked 
for the layers to be considered. 

Then the equation of turbulent diffusion for heat becomes: 



(4D =-v.(7„e + 



V ^ 



Also 









Close to the ground the following approximations are valid: 

-i-ir • . - V'VhO 



1 



Hereinafter the bars over © and T will be omitted but understood. 
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2 i *'0 rr'jp"' 



bno . hnuor(?^iOf>M 



■1 '.. j \iU'A. -'.i vilvi-^D^bnoo v'oba lo Jne;-j;ricoo e.iT 



•\ ■ 

bovi^ - j: .•.■ ■' I .ir.oiiT9v eili ;i; faa;‘ JneludjLiJ snt ci C a''0;'v 
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. t-: !ovii(rvvo> 
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haJi.q <s hnti J- :rs ^'A to sc-blnv '/biiJr. airi': ;il 

'■ , ' 'tC' . 3.'''3 n» eck-v svi-iov/i sit: \'*i r,evio Q rilr-v (I) ,p6 n ■!;■>;. - a 

bssilki; •- : o/Ja ’;* ^nkvciiot sdj ia 'el ;>rTins5j;.' 9 ;lt evods asiilov' lol .cxak.r* 

: . xav vL'ft.iii; oi .70 SK’nariOXS 2T!'=t..idrii f [i5.*noai'(Oii 

-nc ‘ ' '-•■ I ,no,L*oev'''5. J.-..oxtja'v- biw .isoijv&v sdl rtX leriJ oj bB'fajr'Go 

30/57 a. ' '5 c13 J;. 'X 1 -; k-.inoo nod oJ r:o2;7i-iqmoo n* ocf oJ DOiaOi,-. 
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b9>(co' , V 'c -■.. yilcnab io noi^^6h..oY-'Triclsd flcira od.3 bni> , nok'M’.hTno i 

. Laiie'bic.noo sci o+ nioyRi s.i.‘ 'rod 
^d•noua■J. .ToOii jo^i nokidah ine.' oaii/j ic iioil&nqo srfj next'' 



u * ' ' \ 

■ * ' k 









O'.L-. 









: 

) 



:biii5v Gijci anoitGmixoiaqB r’.'Tivv^oIIot erli bnucny oj eeun^ 



J •' V 






'»d V 



. rx/Oj'- ;obni. 'oc bo’JXi'ao ad ilr,v T bno Ta''o -mad odJ Tsliont/j'sH' 



Nov/ raturning to eq. (2) and substituting eq. (3) vve get; 









( 4 ) 



In firne difference form this becomes: 






A 



A-?,.: 






-rcn.-jpo5ina terms and solving for K£ gives: 

(ilLs (« 
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3 « 1 . The water vapor problem . 

The water vapor radiation problem was approached using a method 
prer.ented by D, L. Brooks /!_/ and based upon an adaptation of the 
flux divergence equation of Bruinenberg: 




'^*T 

where q is the specific humidity at the reference level at which 



is to be 



computed; the Elsas ser pressure correction factor; 



■4 



and /i.f'S ^ /is the change in the Stefan-Boltzmann black body flux between 



■successive levels of temperature T in the sounding. The summation 
index n denotes successively higher levels in the atmosphere above the 



reference level and the index m the layers below the reference level. 

The partial derivative ^^is the slbpe of the curve of emissivity of water 

— j- “4 

vapor, plotted against corrected optical depth and the term <?' ■ is 

eV', hunted at the top of the sounding. 

The £actors^i>are evaluated as finite difference slopes from 

dW 

Brooks’ curve^. In addition to the curve Brooks / ly also presents 
t ibiilar values of the mean slope of the emissivity curve for small incre- 
ments of optical path-length . 

If the convention in using eq. (7) is that^i is always positive 

r) 

v«hen temperature decreases with height above the reference level, cooling 



“Not shown here but presented in ^ 1/. 
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at the reference level results from the layer contributions above » 

Sirilarly v/.hen the temperature from the ground decreases up to the ref-= 

erence level, the contributions from below yield warming at the reference 

/ 

le“’el„ The second term on the right, representing radiation divergence 
from the top of the sounding, gives a cooling contribution at the reference 
lerelo The actual conventions used in working with the tables are some- 
vhat different, however, and were designed for maximum efficiency in 
m.*chirie computation. 

In the application of the method, the initial step is the computation 
of the corrected optical depth , w , given by“ 




wh.&re pQ and p. are the layer boundary pressures and p^ is a standard 
pressure, 1000 mb, 

■£q. (3) is approximated for a given layer by using mean values of 
q and p , and a finite pressure increment Ap: 




When increments of optical path length. Aw, are calculated for 
each layer, values of ^^may be obtained by interpolation in the table, 

d vV 

using values w and w + A w in entering the table . 

When this work w'as begun the decision was made to perform the 
computations by hand methods rather than by electronic digital computer. 
This decision was in part based on the fact that use of the then available 
computer, the National Cash Register Company Model 102A, would have 
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necessitated lengthy computation times exclusive of the time required 
for programming. This observation is based on the fact that there are 
over 1000 entries in Brooks' table and only about as many storage cells 
in t> e computer. If the 102A had been used it would have been necessary 
to store the entire table on magnetic tape, which would virtually have 
eliminated the computation time -advantage of the computer. This diffi- 
culty may now have been overcome with the advent of the Control Data 
Corporation model 1604 computer. The larger storage capacity (over 
32,000 cells) and extremely high speed of this new computer would lend 
itself well to Brooks' tabular method. 
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3.2. The ';arbon ■'iioxide * adiation problem . 

For the carbon dioxide part of the radiation problem it was decided 
to employ a simple set of working equations previously used by Brooks 
/i./ in stratospheric investigations of carbon dioxide radiative transfer. 
The downward flux of infrared radiation at a level in the atmosphere 
and for a given frequency band (12.2 - 17.5yw for CO 2 ) is given by the 
Schwarzschild equation integrated over all angles: 

?\/{^\) ~ (O, n/,) 4- )_ Bi A L,-r^ (10) 

C-l 

where /-/•{tv/.)is the downward flux of diffuse radiation at the reference 
level, where ft/- k/', ; w is the optical path length of CO- above the level; 

F I (0) is the flux down from the top of the atmosphere; (0,Wj) is the 
transmissivity of diffuse radiation by a slab of optical path w^; and 

= . 185c 5'T^ is the total black body flux^ radiated by CO 2 at the approp- 

J-T- 

riate mean temperature of the i layer. 

The summation index i runs through n layers from the top of the atmosphere 
to the level in question. is the increment to the flux transmissivity 
at the level w^ contributed by each layer. 

Similarly the upward flux F I' (w^) is: 

2T^ ( F'p- (10) 

where F f (w ) is the upward flux at the ground; and 'TC- (w, ,w_) is the 
o /■ 1 o 

transmissivity (assumed independent of temperature) of the slab of 

^ *A carbon dioxide contribution equal to 18.5% of the black body radiation 
is quoted by various authors . 
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optical thickness w^-wj , or from the ground to the reference level. The 
summation j=l through m runs from the ground to the level wj. 

In using these equations it is assumed that each layer can be app- 
roximated by an isothermal layer whose absolute temperature is the actual 
mean fourth rpot of the temperature of the layer raised to the 4th power. 
Since the soundings do not run to the top of the atmosphere, F (r (0) was 
computed at the top of each sounding. 

The net flux at the reference level (w^) is then simply: 

~ — h 1/ L W,) ( 12 ) 

and the flux divergence, , is taken as being the difference of net 
flux at two adjacent levels , the assumption being that the net flux 
divergence is uniform between the two layers . The hourly-temperature 
change in terms of the flux divergence in the layer is then: 

(•^] = 4,1 K iO AJ. (13) 

2 

where flux is in cal/cm / min and pressure is in millibars. 

For the determination of flux -transmissivity 77= » several choices 
were available as to method. Callender /4_/^ives formulas ior based on 
a breakdown of the 12.5 - 17.5 band into five smaller wave length 
intervals. Burch and Howard /3_/, and Howard, Burch and Williams 
/.B_/ give two empirical formulas based on laboratory data for the total 
absorption of the 12.2 - 17.5 yu band: 

"Strong" fit: -h S'b K/ 4 - 47^ 

"Weak" fit: /Au = 3,/<^ (p^/>y (15) 
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— ^ 4i> di; ^ Al/ 

It follows then that A = ~ ' 

and ^ J — A 

In (14)^ (15) w is absorber concentration in atmosphere-centimeters; p 
is atmospheric pressure in mm Hg; and jS is the partial pressure of CO 2 
(assuming a constant ratio by weight, .00046, of the atmosphere.) 

The Howard et al expressions (14) are for the "strong" fit or rela- 
tively large values of w and (15) for "weak" fit or relatively small values 
or w. The cut-off value given is 20% absorption, the log expression 
applying for 20% and the exponential expression applying for 20%. 

Both of the methods discussed above have the disadvantage of 
extremely lengthy computation-time for the data under consideration unless 
an electronic digital computer is used. Furthermore, the empirical equa- 
tion of Howard et al has the disadvantage of a discontinuity between the 
two equations near absorptivities of 20%. A third and faster choice pre- 
sented itself in the form of a curve giving fractional absorption A vs pres- 
sure X path length (mb-cm-atmos) given by Brooks l_2_J , and credited by 
Brooks to Burch and Howard /_ (see fig. 2.) This curve had the advan- 
tage of affording greater facility in conjunction with hand computations , 
and it had no obvious discontinuities . A difficulty presented itself, 
however, after the CO 2 computations were well underway when it was 
discovered that use of the Brooks curve yielded positive CO 2 radiative 
flux divergences near the ground where large negative water- vapor flux 
divergences were obtained. Since this occurred near the level where 
absorptivities of 20% might be expected, the difficulty was overcome by 
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changing from the curve to the Howard et al "weak" fit expression for 
curve absorptivities less than 20%. 

Although Brooks credits the curve to Burch and Howard / 3_/, it is 
not entirely clear how he obtained his curve. The Burch and Howard curve 
is similar to fhat of Brooks except that their ordinate is percent trans- 
mission and the abscissa parameter is absorber concentration (cm-atmos) . 
The Burch and Howard curve was obtained from data based on the use of 
the two empirical formulas , the discontinuity being smoothed out in the 
construction of the curve. Thus, regardless of how the smoothing was 
done in obtaining Brooks' curve it is clear that it is based on the empirical 
expressions of Burch and Howard. On this basis, changing from the curve, 
fig. 2, to the "weak" fit formula where necessary is justified, the discon- 
tinuity being accepted. 

For the sake of comparison, values of percent absorption due to the 
"weak" fit formula are plotted in fig. 2 along with Brooks' curve. In the 
latter, a value of the "strong" fit formula corresponding to an curve absorp- 
tion of 20% is shown. It appears that the smoothing in the region of 
10-20% absorption favored the "strong" fit curve. 

Since the fractional concentration of CO 2 by weight in the atmos- 
phere is nearly constant with height, the amount near any pressure level 
over a given vertical distance may be expressed in terms of a pressure 
difference. Thus every 4.2 mb of atmospheric air is equivalent to 1 
cm-atmos of carbon dioxide (ie , a layer of CO 2 of thickness 1 cm if 
reduced to S.T.P.). Taking the diffuse radiation through a slab to be 
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2. Percent absorption of terrestrial radiation 
to CO 2 in the 12, 2-17. 5 / 1 band (after Brooks) 








\ 



1.67 times that through a parallel radiating column, the equivalent 



optical path is; 




This is the absorber-concentration in the Burch and Howard expression. 

V. 

For use with Brooks' curve, w is linearly pressure-weighted giving w' as 
follows: 



with w' now in mb-cm-atmos and p the mean pressure of the layer. 
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4. 



Computations . 



4.1. The water vapor radiative temperature change . 

The radiative temperature changes due to water vapor were com- 
puted by use of the method discussed in section 3.1. Required are the 
knowledge of ^temperature , pressure, pressure differences, and relative 
humidity at the given reference levels, the reference levels being taken 
as 8, 17, 35, 100, 165, and 240 meters. Here^ conditions at 16.0 meters 
are assumed to apply at 17 meters , so that in some cases, in order to 
smooth out artificial inversions in the temperature data between 16.0 
and 17 meters, it was necessary to shift the MIT temperatures so that 
the 16.0 meter MIT temperature would correspond to the 17 meter L-20 
temperature . 

The pressure, and layer pressure differences were computed by 
use of the hypsometric equation together with the known temperature 
profile . Values of specific humidity were obtained from the aerograph 
humidity data, radiosonde data, and station humidity. Unfortunately 
instrumental failure left a hole in the micrometeoro logical humidity data 
between the station values and the aerograph soundings . For levels 
between station values and aerograph soundings, linear interpolation 
with height was assumed valid in arriving at values of q for the reference 
levels. As noted, radiosonde data was used to extend the data above 
the maximum level of L-20 flight. This was done by "hooking" on to the 
nearest radiosonde reporting level after the highest reported L-20 data 
were exhausted. 
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The computed water vapor radiative temperature changes are 
given in table 1 . Note the strong radiative warming tendency in the 
lower layers at midday , tending to change to cooling in the neighborhood 
of 100 m. This effect has already been noted by Cooley /5_/. 



Table 1. Infrared radiative temperature change (°C/Hr.) 
due to water vapor (times in CST) . 

Reference Times 



level (meters) 


0635 


08 35 


1035 


1235 


1435 


1635 


8 


-0.13 


1.25 


2.00 


3.05 


2.23 


0.71 


17 


-0.07 


0.12 


0.17 


1.87 


0.92 


0.10 


35 


-0.02 


0.14 


0.20 


0.13 


0.50 


-0.07 


51 


-0.09 


0.09 


0.07 


-0.11 


0.10 


-0.14 


100 


-0.14 


0.03 


0.02 


0.22 


-0.01 


-0.08 


165 


-0.19 


0.09 


-0.03 


0 


-0.01 


-0.01 


240 


-0.17 


-0.20 


-0.05 


-0.03 


-0.06 


-0.05 
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4.2. The radiative temperature change due to carbon dioxide . 

The CO 2 radiative temperature changes were computed by the 
method of section 3.2 and are shown in table 2. Since the mean black 
body flux cT T for layers was needed in the GO 2 method, plots of 

versus height were made for the layers from 0 to 16 meters, and 0.185x 

.-7-4 4 

dT approximated by equal-area averaging . Above 16 meters, 0.185(5' T 

was simply averaged between levels . 

Since the method yields layers between reference levels, 

reference-level values were interpolated. On the average, CO^ radiative 

temperature changes were found to be approximately 23% of the water 

vapor radiative changes, a result which also confirms that of Cooley /_S_/ 

in these layers. Note also the tendency for the radiative temperature change 

to change sign near 240 meters in a manner similar to that of water vapor. 



Table 2» Infrared radiative temperature change 
due to carbon dioxide (°C/Hr.) 



Reference 
level (meters) 


0635 


0835 


Times 

1035 


1235 


1435 


1635 


8 


-0.10 


0.17 


0.26 


0.49 


0.31 


0.13 


17 


-0.02 


0.06 


0.15 


0.26 


0.25 


0.08 


35 


0 


0.03 


0.07 


0.09 


0.10 


0.01 


51 


0 


0.02 


0.05 


0.06 


0.06 


0 


100 


0 


0.02 


0.02 


0.02 


0.02 


0 


165 


0 


0.01 


0.01 


0.01 


0.02 


0.01 


240 


CD 

0 

0 

1 


0 


0.01 


0 


-0.01 


0 
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4.3. The local rate of change of temperature. 



The local rate of change of temperature was obtained from the MIT 
data at 4.0, 8,0, and 16,0 meters and from the aerograph temperature 
data at 35, 51, 100, 165, and 240 meters. Since the observation times 
were centered bi-hourly at 35 minutes after the even hour, ^%t:was taken 
over a two hour period and centered at a time 35 minutes after the odd 
hour. For example, the local change from 0635-0835 CST was centered 
at 0735 CST. The aerograph temperature changes were obtained in a 
like fashion except that some difficulty was encountered in the use of 
this data, A plot of L-20 data with height yields two curves since data 
is taken on both ascent and descent of the aircraft. As an illustration, 
fig, 3 gives the L-20 temperatures for the 1235 sounding. Where data 
for the given reference level were missing smooth-curve interpolation 
was used. This would give an estimated ^ two-hour period either 

from the ascent or descent curves or both. Where doubt existed as to the 
representative values of ^^t^i.e. whether to use the ascent or descent 
curves) , either an average value from both ascent and descent curves , or 
a subjective choice of the most reasonable sounding was used. Since all 
periods were taken as two hour periods , an error in^'^t'^^y possibly 
exist due to a displacement in time away from the supposed observation 
times. This type of error is estimated to be no more than 5% however, and 
should be random in nature. Table 3 presents the observed (or estimated) 
values of ^ . 
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Table 3 . Average local rate of change of temperature (°C/Hr) 



Reference Times 



level (meters) 


0735 


0935 


1135 


1335 


1535 


8 


3.24 


2.27 


1.61 


1.07 


0.12 


17 


2.98 


2.23 


1.55 


1.14 


0.15 


35 


2.08 


2.67 


1.41 


0.80 


0.45 


51 


1.58 


2.68 


1.49 


0.67 


0.35 


100 


0.61 


2.61 


0.87 


0.80 


0.40 


165 


-0.20 


2.15 


0.93 


0.90 


0.37 


240 


0.15 


1.42 


1.05 


0.97 


0.33 
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4.4. Advective temperature changes . 



Generally speaking^in computing the assumption will be made 
that advection is zero. However, it may be shown from the quoted geo- 
strophic winds /lO/ that this assumption is not quite valid for the period 
in question. Haltiner and Martin / 7/ have shown that the layer mean 
geostrophic advection is given by: 

- V 5 - Vp T, = if l\l\ JU, ^ (18) 

where is the layer mean geostrophic wind; Tv is the layer mean virtual 
temperature; and /I A is the cross-hatched area of the geostrophic wind 
hodograph in fig . 4 . 

Now computing the geostrophic mean advection for the layer 950-900 
mb for 0900 CST on 7 September we get 0.24°C/Hr. Since the actual wind 
component perpendicular to the mean isotherms is about 17% of the mean 
geostrophic, the advection term of 0.24°c/Hr. will be reduced to 
0.041°C/Hr. Since the wind speeds increase in the afternoon it is pro- 
bable that a somewhat larger advection will take place at this time. Un- 
fortunately no other useful geostrophic components are given in the data 
other than for 0900. This consideration of advection, therefore, requiring 
more extensive investigation will not be pursued further in this paper^ . 

^'Except for the 0935 sounding where the advection computed for 0900 
will be used in the determination of a Kj^^for comparison with a 
computed without the advection consideration. 
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4.5. The coefficient of eddy conductivity. 

Having computed /4T) it will be necessary to obtain /^4 ya-) 

Values of /AT. x V>) were obtained by adding 0.0098°C/m to A“r 

U3 ' ITi 

To obtain ^ above 17 meters , plots of aerograph temperature with 
A3 

height were m^de on semi-log paper (see fig. 3). 

Similarly for values up to 17 meters^ plots of MIT micrometerological 

temperatures with heights were made on linear graph paper and a smooth 

curve drawn through the points . Values of AT were then obtained from 

a2 

the applicable curve by measuring the slope of the chord centered on the 
level in question. For vertical gradients at 17 meters^ and below, one 
meter intervals were taken. For the levels 35 meters to 165 meters, IO- 
meter intervals were used, while 20-meter intervals were used above 
165 meters . 

Returning to eq. (6) and omitting advection, one obtains: 

[W - ^ , 20 ) 



k,- 






Because of the fact that was centered after the odd hours and 
the radiative temperature change was centered after the even hours, it 
was decided to compute a K^j centered after the odd hour by doing the 
following: 

(1) Averaging /AT \ for two adjacent after "the even hour" 

soundings., giving a fXT) centered after "the odd hour". 

\A±:/RAJ) 



^'The 16 meter change was assumed to apply at 17 meters. 
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(2) Repeat step (1) again for 4J~ +- ^ to obtair 




values of which are listed in table 4. 



The previous steps yield the final recursion formula 






( 21 ) 



where , K 2 are also centered after "the odd hours" . 

Table 4. Average lapse rate of potential temperature, (°C/meter). 



Reference 
levels (meters) 


0735 


Time 

0935 


1135 


1335 


1535 


4 


0.0088 


t— 1 

r— 1 

0 

t— 1 

0 

1 


-0.1301 


-0.1292 


-0.0992 


8 


0.0173 


-0.0277 


-0.0402 


-0.0482 


-0.0287 


17 


0.0218 


-0.0132 


-0.0252 


-0.0202 


-0.0082 


35 


0.0138 


-0.0137 


-0.0047 


0.0072 


-0.0082 


51 


0.0183 


-0.0057 


-0.0042 


-0.0052 


-0.0062 


100 


0.0178 


0.0023 


-0.0042 


-0.0032 


-0.0042 


165 


0.0208 


0.0073 


-0.0027 


0.0008 


0.0003 


240 


0.0270 


0.0143 


0.0014 


0.0003 


0.0018 



As an illustration of the form of this equation for obtaining at 17 



meters, the equation is; 




with At 2 = 17 meters - 4 meters = 13 meters. 
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Now a starting value of K at four meters , with temperature changes 

H 

at eight meters, yields a value of Kj^ at 17 meters. Next to obtain the 35- 
meter Ktr a starting value of K„ at eight meters with temperature changes 

IT. 

at 17 meters, is required. After this the computed value of Kj^ at 17 meters 
is used to get a value at 51 meter s^ and the scheme is continued in leapfrog 
fashion until all values have been obtained. 

As previously noted in section 1, table 5. a of /lO/ gives values of 
Q (meal cm ^ min ■'■ ) . These values were obtained from the correlation of 
rapid fluctuations of temperature and vertical motion using one-second 
samplings of fast response probes together with the Reynold's formula: 

Q = 

As a general principle, following Priestly /H/, the turbulent flux 
of heat may be considered constant with elevation in shallow layers above 
about one meter. The reported values, table 5. a of /l^/, seem to sub- 
stantiate this result up to 12 meters. Now values of Q given in table 5.2a 
of ^0_/ are for 3.0, 6.0 and 12.0 meters. Averaging these to get represen- 
tative values of Q in the layer 3-12 meters gives starting values of Kpj 
and Kpj g obtained by using eq. (1). 

One further complication exists, however, in the 1235 and 1435 
soundings for which no MIT values of Q are given. Here it was decided 
to estimate values of Q from a comparison of the MIT values of Q and 
theoretical model values of Qq at the earth-air interface. Two such 
theoretical models^ were available with values for the 1235 and 1435 

1 . - - 

For a discussion of these theoretical models see Chapter 7 . 3 of 
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■ n^- 







. c i- . t s ; one due to Suomi AQ_/ and one due to Lettau 

Considering the entire seventh observation period^Lettau's values 
of Qq were found to correlate more highly with the reported values of 
turbulent heat transfer Q than did the values of Suomi. On this basis, 
a ratio of the.average observed Q to the Lettau theoretical value Qq for 
the 1035, 1235, and 1435 observations for aU observation periods at 
O'Neill, Nebraska was calculated. 1035 values were included to give 
enough values on which to base this ratio. The ratio was then applied 
to the Qq of the Lettau theory for 1235 and 1435 so that and Q 0 for 
these times could be obtained. Computed values of Kj^ are given in 
table 5 . 
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Table 5. Computed values of K„, (cm^/sec.) 
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5. 



Results . 



5.1. Kj^ as a function of height . 

Considering the values of given in table 5 up to the point where 

there is a sharp drop off in values with respect to height^ , the computed 

values of «how a strong linear trend when plotted on log-log paper. 

To substantiate this linear trendy correlation coefficients of log 

Kj^ and log "Z were computed and regression lines of best fit were obtained 

and plotted for the data for each sounding (see figs. 5-10). 

Now a straight line plot of K„ with Z on log-log paper yields K„ 

n hi 

as a function of exponential height according to the general equation: 




Computed equations for K„ as an exponential function of height are 

ti 

presented in table 6. 

Table 6. Kj^ as an exponential function of height. 

Range of Heights on 



Time (CST) 


Equation 


which the equation is based 
(meters) 


0735 


Kn = 2: ' 


4 - 320 


0935 


2- ' 


4 - 320 


0935^ 




4 - 320 


1135 


-2Cc( H 


4 - 165 


1335 


2. 


4 - 165 


1535 


/Ch - /5'00 2L 


4 - 100 


Negative values, and values above 165 meters for the 1135 sounding, 
were excluded . 

’Advection neglected in the computations except for the second result at 
0935, for which the 0900 computed advection was assumed to apply at 0935. 
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Fig. 5. Ky as a function of height for the 0735 sounding. 
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fuiictiun of height for the 0935 ? sounding (-without advection.) 
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a function of height for the 0935 sounding (with adveotion.) 





Fig. 8. Kn as a function of height for the 1135 soxmding. 
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as a function of height for the 1335 sounding. 
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Fip, 10. Ky as a function of height for the 1535 sounding 




It IS to be noted that where advection was considered (0935 case) 

the result was not greatly different from that in which this parameter was 

neglected. The difference, however, indicates that the Kj^- values, 

with advection included, were slightly greater in the lower portion of the 

layer, and smaller in the upper portion. 

It should also be pointed out that in most cases, the values of K„ 

rl 

computed for elevations above 165 meters showed a decrease below the 

value at the previous level. In several cases, the values of K„ were 

hi 



negative. This was most likely due to observational errors in , 

when was very close to the adiabatic lapse, a condition which 

occurred in the neighborhood of 165 meters for the midday and afternoon 



runs . 
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5.2. K„ as a function of stability. 
ri 

It would have been appropriate to compare values of and m 
to a stability parameter such as the Richardson number. Since appropriate 

wind data for such a comparison were not available, it was more feasible 

_ 2 
to correlate K and m with a significant part of the numerator of the in- 
H 

tegrated Richardson number^ A 0 . This was done and the results are tabu- 
lated below. 



Table 7. Linear correlation coefficients of A0with K andA0 with m 



at indicated heights 



H 



K (meters) a 0 computed Correlation 

^ over layers (meters) coefficients of 

Kh with A0 



Correlation 
coefficients of 
m with A0 



4 


2-4 


0.26 


CO 

o 


8. 


2-8 


-0.77 


-0.71 


17 


2-17 


-0.96 


00 

(X) 

o 

1 


35 


2-35 


00 

CO 

o 

1 


1 

o 

00 


51 


2-51 


-0.76 


-0.75 


100 


2-100 


-0.67 


-0.75 


Considering values of A© versus 


Kpj at heights from four to 100 



meters for all soundings (sample size 30) a correlation coefficient of 
-0.43 is obtained. <• 



1 

2 



Values obtained from equations given in table 6. 
Here A© = ©(z)- 0(:^) 

where ^ is variable and ©(2) is the value of 0 at Zi 



= 2 meters . 
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6 . 



Conclusions . 



Although negative values of are somewhat disconcerting, other 
investigators have been confronted with this anomaly at much lower levels 

As mentioned earlier, factors possibly contributing to these negative 
values are chiefly non-representative lapse rates due to errors in the 
measurement of temperature, especially as the lapse rate of temperature 
approaches the dry adiabatic. Still another possibility arises when one 
considers the theory of Priestly and Swinbank concerning forced versus 
free convection (see for example p. 261 of /j _/) . 

From an inspection of the computed power-law profiles for listed 
in table 6, it appears that m increases inversely with stability, the maxi- 
mum value of m occuring in the early afternoon. Although it is doubtful 
whether the correlations listed in table 6 can be proved statistically sig- 
nificant on the basis of such a small sample size (five) , nevertheless it 
is felt that a strong physical significance is inherent in the repetative 

nature of the results, and that K has a strong inverse dependence on 

H 

stability. This statement is given support by the results of a correlation 
of the A0 parameter with all values of Kj^ represented in table 7 (sample 
size 30) where a statistically significant^ correlation coefficient of -0.43 
is obtained. Further it is apparent, from an inspection of table 7 that this 
correlation (-0.43) could be further improved by the elimination from con- 
sideration of values of at four meters where a positive correlation of 

with A0 is obtained. Considering the positive correlation of K„ 
n , i H , 

^ ‘At the five percent level of belief. 
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with A© and the negative correlation of K„ n with A© , it seems likely 
that the method of computing at four and eight meters may have affected 
the overall correlation. It must be remembered that the assumption made 
in this connection was that is independent of elevation from one 

to 12 meters. At higher levels, the computational procedure gave con- 

a. 

siderably larger Kj^-values which then varied consistently (if not '’signi- 
ficantly") with stability. 

In summary, daytime-values of obtained under the conditions of 
the Great Plains Turbulence Field Program /lO/ , tend to increase with in- 
creasing lapse rate, and with increasing height up to at least 100 meters. 
Above some such level Kjj seems to decrease. It is possible that a 
specific function may exist, but the data was not sufficiently 

detailed to establish such a result at this time. 
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